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Abstract. We reexamine arguments advanced by Hayashi & Matsuda 
(2001), who claim that several simple, physically motivated deriva- 
tions based on mean free path theory for calculating the viscous torque 
in a quasi-Keplerian accretion disk yield results that are inconsistent 
with the generally accepted model. If correct, the ideas proposed by 
Hayashi & Matsuda would radically alter our understanding of the na- 
ture of the angular momentum transport in the disk, which is a central 
feature of accretion disk theory. However, in this paper we point out 
several fallacies in their arguments and show that there indeed exists a 
simple derivation based on mean free path theory that yields an expres- 
sion for the viscous torque that is proportional to the radial derivative 
of the angular velocity in the accretion disk, as expected. The deriva- 
tion is based on the analysis of the epicyclic motion of gas parcels in 
adjacent eddies in the disk. 
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1. Introduction 



Accretion disks around a central compact object are ubiquitous in astrophysics. In 
such systems, the plasma accreting onto the central object typically takes the form 
of a swirling disk. The disk-like nature of accretion is due to the angular momen- 
tum originally carried by the accreting material, and also to the spin of the central 
object, which provides a preferred axis of rotation. At a given radius in a quasi- 
Keplerian accretion disk, the average motion of the gas parcels is an approximately 
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circular Keplerian orbits, combined with a radial sinking motion as the parcel grad- 
ually moves deeper into the potential well of the central mass. The radial velocity 
that is far smaller than the Keplerian (azimuthal) velocity. A detailed description 
of quasi-Keplerian accretion disks can be found in Pringle (1981), among several 
other references. 

The specific angular momentum (angular momentum per unit mass) carried by 
a parcel of gas in a quasi-Keplerian accretion disk is very close to the Keplerian 
value \J GMR, where R is the radius measured from the central object with mass 
M. It follows that a parcel of plasma has to lose angular momentum in order to 
sink from a larger orbit into a smaller one and eventually cross the event horizon. 
In the absence of strong winds or jets, the angular momentum must therefore flow 
outwards in such an accretion disk in order to enable accretion to proceed. In 
disks of turbulent fluid, the flow of angular momentum is due to the exchange of 
plasma parcels in neighboring annuli with different angular velocities. This fluid 
picture may describe certain astrophysical cases such as the rings of Saturn, but in 
many other cases of astrophysical interest, the coupling between adjacent annuli is 
provided not by "parcel interchange," but rather by the magnetic field. In either the 
fluid or magnetohydrodynamical (MHD) scenarios, the angular momentum flow is 
expressed as a torque between neighboring annuli. Our focus here is on the fluid 
picture, and on the validity of the various attempts to analyze quantitatively the 
angular momentum transport associated with the interchange of parcels. 

There are several previous papers in which derivations of the viscous torque in 
fluid dynamical situations are developed based on simple, physically motivated ar- 
guments (e.g., Frank, King, & Raine 1985, 1992, 2002; Hartmann 1998). However, 
some of these do not yield an expression for the viscous torque that is proportional 
to the radial derivative of the angular velocity. This is a problem because the vis- 
cosity is fundamentally due to the "rubbing" of matter in adjacent radial annuli in 
the disk, and consequently the viscosity should vanish in the case of solid body 
rotation with Q(R) = constant, where Q(R) is the angular velocity in the disk 
at radius R. Hayashi & Matsuda (2001; hereafter HM) recognized this point, and 
attempted to clear up some of the confusion by carefully examining the previously 
published derivations of the viscosity that were based on the mean free path ap- 
proach. They concluded that the mean free path approach inevitably leads to an 
inward rather than outward flow of angular momentum in the disk, which is un- 
physical. However, we argue that the reasoning of HM was flawed because they 
did not consider the epicyclic nature of the parcel trajectories. We present a simple 
derivation based on mean free path theory, combined with the actual epicyclic mo- 
tion of the gas parcels, that in fact yields a physically reasonable expression for the 
viscous torque between neighboring annuli in a quasi-Keplerian accretion disk. 

We discuss the standard fluid dynamical formulation of viscous torques in an 
accretion disk in § 2. We then discuss the problems with the various previous 
derivations that attempted to utilize a mean free path approach to compute the vis- 
cous torque in § 3. In § 4 we present a simple physical derivation based on analysis 
of the epicyclic (ballistic) motion of gas parcels in adjacent eddies, and we demon- 
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strate that this approach yields the expected form for the torque in terms of the 
gradient of the angular velocity. We present our final conclusions in § 5. 



2. Standard fluid dynamics treatment of viscous torque 

The equation of motion for a viscous, incompressible fluid can be written as (Lan- 
dau & Lifshitz 1987) 

where the force density F is defined via the viscous stress tensor as follows: 

p _ daik 

( dvi dv k \ 
Oik = ~P Oik + r}[ j: h ^ — I . (2) 

The fluid velocity is denoted by v and its density and pressure by p and p, respec- 
tively. The coefficient of dynamic viscosity (g cm -1 s _1 in cgs units) is denoted by 
7] and we have used the Einstein summation convention in equation (2). For the spe- 
cial case of a thin, azimuthally symmetric accretion disk, the only non-negligible 
component of the viscous stress tensor in cylindrical coordinates (R, <j>, z) is 

OR4> = -VR-Qft- (3) 

The viscous stress (force per unit area) is thus directly proportional to the radial 
derivative of the angular velocity. This is by far the cleanest and most rigorous way 
to derive the azimuthal equation of motion for a quasi-Keplerian accretion disk 
(see, for example, chapter 1 of Subramanian 1997). It is implicity assumed that 
the dynamic viscosity i] arises out of local effects; i.e., due to momentum exchange 
between neighboring annuli of the accretion disk, as with molecular viscosity. It is 
well known (Shakura & Sunyaev 1973; Pringle 1981) that molecular viscosity is far 
too small to account for angular momentum transport in accretion disks around ac- 
tive galactic nuclei. Identifying suitable candidates for the microphysical viscosity 
mechanism operative in such disks is a subject of intensive research. 



3. Simplified treatments of viscous torque in accretion disks 

Pringle (1981) derived the azimuthal component of the equation of motion in ac- 
cretion disks starting from first principles in a simple, physically motivated manner. 
The equation reads as follows: 

d (S R 3 n) d(Zv R R 3 n) \_dQ_ 

dt + d~R ~~27aR' () 
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where £ is the surface density of plasma in the disk, v R < is the radial accre- 
tion velocity, and Q > is the torque exerted by the material inside radius R on 
the material outside that radius. Throughout the remainder of the paper, we shall 
assume that the disk has a steady-state (time-independent) stucture, although this 
is not essential for our results. In the standard approach introduced by Shakura & 
Sunyaev (1973) and adopted by Pringle (1981), the torque Q is related to Q via 

Q = 47rR 2 Ha R4> = -2vri? 3 S ^ , (5) 

where H(R) is the half-thickness of the disk at radius R and v = rj/p is the 
kinematic viscosity coefficient. The stress and torque are therefore proportional 
to dVL/dR, in agreement with equation (3). This prescription has been applied in 
many disk stucture calculations. In particular, it has been shown recently (Becker 
& Le 2003) that fully relativistic and self-consistent models for hot, advection- 
dominated accretion disks can be constructed by applying the Shakura-Sunyaev 
viscosity prescription throughout the entire disk, including the region close to the 
event horizon. 

A number of authors have attempted to confirm the general form of equation (5) 
by using simple physical arguments. However, several of these derivations have er- 
rors in them, and they do not always result in an expression for Q that is proportional 
to dft j dR, as pointed out by HM. We briefly review the relevant derivations below, 
and we also point out errors in the approach adopted by HM. We then present a 
new, heuristic derivation of equation (5) that is based on a careful analysis of the 
ballistic motion of two parcels as they exchange radii. This derivation leads to the 
expected conclusion that Q oc d£l/ dR. 

Although the argument given by HM is rather indirect, their main point can 
be understood through a simple examination of the angular momentum transport 
resulting from the interchange of fluid elements in a disk that is rotating as a solid 
body, i.e., with = S7o = constant. In this case, the angular momentum per 

unit mass, denoted by J = R 2 Q(R), is given by J = R? Sl > an d this quantity 
increases rather strongly as a function of the radius R. Hence if two parcels of fluid 
on opposite sides of radius R were exchanged due to some turbulent or convective 
process, then clearly angular momentum would be transported in the inward direc- 
tion, since the blob that was originally outside the annulus will have more angular 
momentum than the interior blob. However, this result is unphysical, because in the 
case of solid body rotation, there is no "rubbing" between adjacent fluid annuli, and 
therefore there should be no torque and no angular momentum transport. Any suc- 
cessful microphysical model for the angular momentum transport in the disk based 
on mean free path theory must somehow resolve this apparent paradox. In the fol- 
lowing sections, we provide a detailed consideration of the reasoning employed 
in the previous published derivations, including that of HM, and we conclude that 
when properly carried out, the mean free path approach can yield a result for the 
viscous torque that correctly vanishes in the case of solid body rotation. 
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3.1 Derivation ofQ by Hartmann (1998) 

The specific angular momentum J m carried by material originating at a radius R — 
A/2 in a quasi-Keplerian accretion disk is given by Hartmann (1998) as J in = 
(R — A/2) — A/2) , where A is the mean free path over which parcels of 
plasma exchange angular momentum. As pointed out by HM, this expression is 
incorrect, and the correct expression should read as follows: 

Jin= { R ~^} 2 Q { R ~l) ■ (6) 

One can write an analogous expression for J out , the specific angular momentum 
carried by material originating a a radius R + A/2, by reversing the sign of A. As 
shown by HM, if one expands Q(R — A/2) to first order in A as Q(R — A/2) ~ 
fi(i?) — (X/2)(d£l/dR), the result obtained for the difference between the specific 
angular momenta is 

Ji„ - Jout = - A ^ {R 2 n) . (7) 

The characteristic time for the interchange of the matter between the two radii is 
At = X/w, where w is the turbulent velocity of the fluid parcels. The total mass 
of fluid involved in the interchange is AM = 2-kRXT,, and it follows that the net 
rate of flow of angular momentum from the inner ring at R — A/2 towards the outer 
ring at R + A/2 (or equivalently, the viscous torque exerted by the ring at R — A/2 
on the ring at R + A/2) is 

AM 

Q = (Jn - Jout) = 2nR S W (J in - J out ) , (8) 

or 

G = -2nRi:[3v^-(R 2 n), (9) 
dR 

where we have set v = w\/ [3, with [5 denoting a constant of order unity. The 
viscous torque in a quasi-Keplerian accretion disk should tend to smooth out gra- 
dients in the angular velocity so as to attain solid body rotation (dfl/dR = 0). 
However, the expression for the viscous torque in equation (9) is such that it tends 
to attain a flow with constant angular momentum, i.e., the "equilibrium" condition 
is d(R 2 Vi)/dR = 0. As HM point out, this expression is therefore unphysical. 



3.2 Derivation of Q by Frank, King, & Raine (1992) 

We next turn our attention to another derivation of the viscous torque given by 
Frank, King, & Raine (1992). At the heart of their derivation is the claim that the 
linear velocity of material at radius R — A/2 as seen by an observer situated at 
radius R is given by 

v ie i(R-^j = (R-£)n(R-£)+Sl(R)±. (10) 
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Based on this, one can write an expression for L m , the rate at which angular mo- 
mentum crosses an annulus at radius R in the direction of increasing R, as 

L in = 2irRVw(R- ^) v rcl (R-^ 



2irRZw 



- 2irR £ w I R 



RSl(R) - R 



A dfl 
2 ZR 



(11) 



where we have expanded 0(12 — A/2) to first order in A to arrive at the final ex- 
pression. We can write an analogous expression for L out , the rate at which angular 
momentum crosses radius R in the inward direction, by reversing the sign of A. To 
first order in A, the torque exerted on the plasma outside radius R by the plasma 
inside that radius is then 



G = L in - Lout = ~2ttR 2 pvX^-(Rn). 

(IK 



(12) 



Since the right-hand side is proportional to d(RQ)/dR, in this case the torque 
will lead to a uniform linear velocity, and not to a uniform angular velocity as we 
require on physical grounds. As HM point out, this expression for Q is therefore 
also incorrect. 



3.3 Correction proposed by HM 

HM claim that this is because Frank, King, & Raine (1992) have used an incorrect 
expression for v re \ (i.e., eq. [10]). They assert that the linear velocity of the plasma 
at (R — A/2) as viewed by an observer at radius R should instead be given by 

v rel (R-^j = (R-£) n ( R -^) -RV(R) + n(R)^. (13) 

By employing this expression for v ie \ and following the same procedure used to 
obtain equation (12), they find that to first order in A, 

g = L in -L out = -27r J R 3 /3^S^. (14) 

dR 

This expression for Q does indeed have the correct dependence on dQ/dR, but 
nonetheless we claim that equation (13) for the relative velocity used by HM is 
incorrect. The correct expressions for the relative velocities are in fact (see, e.g., 
Mihalas & Binney 1981) 

^rel (R ~ ^) = (R ~ ^) « (R ~ ^) - R*l(R) 

v Tcl (R+^j = (r+£) Sl(R+£)-Rn(R), (15) 
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where v Te \(R — A/2) denotes the velocity of a plasma parcel at R — A/2 as seen 
by an observer at R, and and v ie i(R + A/2) denotes the velocity of a plasma parcel 
at R + A/2 as seen by an observer at R. In a quasi-Keplerian accretion disk, 
SI oc i?~ 3 / 2 , and therefore the first velocity is positive and the second is negative. 
These expressions for the relative velocities assume that the plasma parcels at R, 
R — A/2, and R + A/2 all lie on the same radial line through the central object. 
If this is not true, then the expressions will contain additional terms, as Mihalas 
& Binney (1981) show in the context of the relative velocity between the sun and 
stars in our galaxy. Since we are dealing with material transport over lengths of the 
order of A that are very small in comparison with R, this assumption is quite valid 
in our accretion disk application. Using equation (15) for the relative velocities, we 
now obtain 



L in = 2ttR E w [R-- 



L out = 2vri? Sw ( R 



R 



Rn(R) 



~ -2irR 2 Y,w 



A 



9 A 
2 



(16) 



where the final expressions for Lj n and L out are correct to first order in A. We thus 
obtain for the viscous torque 



G = L in - L out = - 2nR 2 £ (3 v^E ( Rn ) ■ 



(17) 



This result is identical to equation (12), and therefore it too is incorrect. 



4. Derivation based on epicyclic parcel motion 

We have demonstrated in § 3 that due to various errors, the previous derivations of 
the viscous torque based on mean free path theory do not yield results consistent 
with the classical theory of viscous transport. We shall now focus on a heuristic 
analysis of the angular momentum transport that combines the mean free path ap- 
proach with a proper treatment of the epicyclic motion of parcels of gas in adjacent 
eddies in the accretion disk. The physical picture is presented in Fig. 1. Let us 
suppose that A represents the mean distance a parcel travels freely before having 
its motion disrupted by interaction with the surrounding material in the disk. The 
disruption in this case involves hydrodynamical interaction, and therefore A is the 
damping length for the turbulence as well as the mean free path for the parcel mo- 
tion. It follows that the turbulence in the disk can be characterized by whirls with a 
mean radius A. 

Next we focus on the interchange of two gas parcels of unit mass, A and B, 
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Figure 1. Angular momentum transport in the disk involves the interchange of two parcels, 
A and B, initally located at radii R+X/2 and R— A/2, respectively, around central mass M. 
They each participate in ballistic, epicyclic motion in eddies of radius A, with A moving 
inward and B moving outward. See the discussion in the text. 



which are initially located at radii R + X/2 and R — A/2, respectively (see Fig. 1). 
The two parcels each "ride" turbulent eddies with radius A, and they start out at 
the outer (parcel A) or inner (parcel B) edge of their respective eddies. Hence 
each begins its motion with zero radial velocity, at a turning point in its orbit. 
Parcel A moves in the inward direction, and therefore its angular momentum is 
sub-Keplerian compared with the average disk at radius R + X/2. Similarly, parcel 
B is super-Keplerian compared with the mean disk at radius R — A/2 and therefore 
it moves in the outward direction. The motion of the parcels is epicyclic as viewed 
from the reference frame of an observer who remains at the starting radius and 
travels along with the Keplerian angular velocity. Since A is the mean damping 
length for the turbulence, each of the parcels will experience ballistic motion over 
a radial length scale comparable to A, after which the ballistic motion (and the 
parcel itself) will be absorbed by the surrounding gas. It follows that, on average, 
parcel A will deposit its angular momentum (J ou t) at radius R — A/2, and parcel B 
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will deposit its angular momentum (Ji n ) at radius R + A/2. Hence the two parcels 
exchange radii during the process. We shall proceed to compute the net angular 
momentum exchange, J; n — J out , by considering the motion of these parcels in 
detail. 

Gas parcel A participates in the motion of turbulent eddy A, which has its 
center located at radius R — A/2 and possesses turning points at radii R + X/2 and 
R — 3A/2. Likewise, parcel B experiences the motion of eddy B, with its center 
located at radius R + X/2, and with turning points a radii R + 3A/2 and R — A/2. 
The specific energy E and the specific angular momentum J of each parcel is 
conserved during the ballistic phase of its motion, until it travels a mean distance 
A and is damped by the surrounding material. For a particle moving ballistically in 
the Newtonian gravitational field of a central mass M, the quantities E and J are 
related by the classical energy equation 

„ 1 2 1 J 2 GM 

E= 2'i + 2W~nf <18) 

where v R is the radial component of the velocity. Turning points in the radial 
motion occur where v R vanishes, so that we can write 

1J 2 GM _1 J 2 GM 

2R\ R l 2R\ R 2 ' U } 

where R\ and R 2 denote the two turning point radii. This equation can be easily 
solved for the angular momentum J as a function of R\ and R 2 . The result obtained 
is 

'2GMR l R 2 \ 1/2 



J = (n^f) ■ <20) 

We can use equation (20) to conclude that the angular momentum of parcel A is 
equal to 



-1/2 



Likewise, the angular momentum of parcel B is given by 



1 A 



3A\ 1 / 2 / 1A 2 V 1/2 



^^l 1 -^^) [ l -4l?) • (22) 

In the spirit of the mean free path approach, we are interested in computing the 
value of the net angular momentum transport, J m — J out , to first order in the small 
parameter X/R. The corresponding results obtained for J m and J out are 



y ()lll = V GMR (l---)+0 



1 A 



R 2 



(23) 
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J m = Vgmr (i + +0 



R 2 



(24) 



The angular velocity Q(R) in a quasi-Keplerian accretion disk is very close to 
the Keplerian value, and therefore we can write 



n(R) 



'GM 



It follows that 



r 2 n [R + 



Vgmr ( i + 



1 A\-V 2 
6R 



(25) 



(26) 



or, to first order in X/R, 

R 2 n 

This also implies that 













4 RJ 



R 2 



(27) 



fi I R - 



R 2 



(28) 



Comparing equations (23) and (24) with equations (27) and (28), we find that to 
first order in X/R, the net angular momentum transfer is given by 



Ju 



X 



6 



J out = nr (l [R - - - R z n (R + - ~ -- R z — 



X 



6 



A 



dn 



dR 



By following the same steps leading to equation (8), we now obtain 

g = 2irRZw (J in - J out ) ~ -2vri? 3 v S ^ , 

Ctrl 



(29) 



(30) 



where the final expression is correct to first order in A and we have set v = wXfi 
so that (3 = 3. Note that this result agrees very well with the Shakura-Sunyaev 
form (eq. [5]). Hence we have demonstated using a simple heuristic derivation that 
the viscous torque is indeed proportional to the gradient of the angular velocity in 
an accretion disk within the context of a mean free path, parcel-exchange picture. 



5. Conclusion 



The derivation presented in § 4 clearly employs mean free path theory, since we as- 
sumed that the fluid parcels travel an average radial distance A before being damped 
by the surrounding gas. We have thus shown that there does exist a simple deriva- 
tion, based on mean free path theory, that yields an expression for the viscous 
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torque Q (eq. [30]) that is directly proportional to the radial derivative of the angu- 
lar velocity, d£l/dR, in agreement with our physical expectation. This expression 
for Q can be used in equation (4) to proceed further in deriving the structure of the 
accretion disk, as in Pringle (1981). Our results provide a simple but important 
unification of the "parcel interchange" viscosity model with the MHD viscosity 
model, in which the coupling that transports the angular momentum is provided 
by the magnetic field rather than by fluid turbulence (Frank, King, & Raine 2002). 
In the MHD model, the torque is found to be proportional to the gradient of the 
angular velocity, dfl / dR, in agreement with the results we have obtained here by 
applying the mean free path approach to the case of fluid turbulence. 
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